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Abstract 

The contribution of the fourth generation quarks to the rare Bs — > and Bs — > 
i^i~'y decays are studied. Constraints on the extended Cabibbo-Kobayashi-Maskawa 
matrix elements |VtbV^*| and |Vt'bV^*g| are obtained from the existing experimental 
result for the branching ratio of the B Xg'j decay. The branching ratios of the 
above-mentioned decays, as well as the photon spectrum in the Bg i^£^^ decay, 
are studied. It is found that the results are sensitive to the existence of the fourth 
generation. 



PACS number(s): 12.60.-i, 13.20.-v, 13.20.He 



*e-niail: taliev@metu.edu.tr 
^e-mail: altugoz@metu.edu.tr 
■l-e-mail: savci@metu.edu.tr 



1 Introduction 



The working and forthcoming high-statistics B-physics experiments at BaBar, BELLE, 
HERA-B, the TeVatron and LHC-B [1] will probe the flavor sector of the Standard Model 
(SM) with high precision. These experiments may provide hints for the physics beyond the 
SM. Already, there have been a large amount of theoretical investigations within various 
scenarios of new physics to calculate the physical observables that will be measured. 

The flavor-changing neutral current process Bs(d) —>■ has attracted a lot of attention 
due to its sensitivity both to the gauge structure of the SM and to the possible new physics 
beyond SM [2]-[9]. Note that the experimental bound on the branching ratio is [10] 

B{Bd i^^i^-) < 6.8 X 10"^ {CL = 90%) , 
B{Bs n^ii-) < 2.0 X 10-*^ (CL = 90%) . 

In other words, measuring the branching ratio Bs(^d) ~^ i~^£~ can impose the stringiest 
constraints on the possible parameter space of new physics. The Bg i^i^'-f decay is 
also very sensitive to the existence of the new physics beyond SM. In spite of the fact 
that this decay has an extra a factor, it was shown in [11, 12] that the radiative leptonic 
Bg — > t^l~^ {£ — e, jj) decays have larger decay rates compared to that purely leptonic 
ones, while for the r Icpton mode widths for both decays are comparable to each other. 

As has already been noted, the flavor-changing neutral current processes Bg — ^ 
and Bg — £+£^7 are very sensitive to the new physics beyond the SM. One of the most 
straightforward and economic extension of the SM is adding the fourth generation to the 
fermionic sector, similar to the three generation case. 

Effects of the extra generation to the electrowcak radiative correction and its influence 
on the low energy physics processes were investigated in many works (see for example [13] 
and references therein). 

In this work we study the effects of fourth generation in the Bg — > t^£~ and Bg — > t^£~^ 
decays. The paper is organized as follows. In section 2 we present theoretical description of 
the Bg — > £'^l~ and Bg — > £+£"7 decays for the fourth generation case. Section 3 is devoted 
to the numerical analysis and conclusion. 

2 Theoretical background 

In this section we will present the necessary theoretical expressions for the Bg — ^ and 
Bg — > £+£~7 decays. 

Let us flrst consider Bg — > decay. At quark level, this decay is described by the 
h s£~^£~ transition for which the effective Hamiltonian can be written as 

-Heff = E Q{n)0,{n) , (1) 

where the full set of operators in Oi(/i) and the corresponding expressions for the Wilson 
coefficients in the SM are given in [14, 15]. As has been noted earlier, in the model we 
consider in the present work, the fourth generation is introduced in the same way as the 
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first three generations are introduced in the SM. The fourth generation changes only values 
of the Wilson coefficients Cr{iJ,), Cq{ij,) and Cio(/i) via the virtual exchange of the fourth 
generation up quark t', i.e., 

Vtb Vts 
Vtb Vts 

where the last terms in these expressions correspond to the t' quark contributions to the 
Wilson coefficients and Vfb and Vfs are the elements of the 4x4 Cabibbo-Kobayashi- 
Maskawa (CKM) matrix. The explicit forms of the C'f^"' can easily be obtained from the 
corresponding Wilson coefficient expressions in SM by simply substituting rrit rrit'. The 
effective Hamiltonian (1) leads to the following matrix element for the b s£^£~ transition 



Ga 



(3) 



where = {pi +^2)^ and pi and p2 are the final leptons four-momenta. It follows from Eq. 
(3) that in order to calculate the matrix element for the Bg t^l" decay, M. is sandwiched 
between vacuum and Bg meson states, from which the matrix elements (0 |s7/i(l — 75)^! Bg) 
and (0 \s(Tni,q"{l + 75)6! Bg) need to be calculated. These matrix elements are equal to 

(0|57'^(1-75)&|5.) = -ifsyB, 
{Q\sa^,q^{l + -i^)h\Bg) = 0. (4) 

Using these results, we get for the matrix element for Bg —>■ decay 

Ga 

M = t-j^Vt,V;gmjBMoimb) V • (5) 

Obviously, we see that the vector leptonic operator l7^£ does not contribute to Bg — > i'^i" 
decay because it gives zero when contracted with Bg meson momentum. 

Let us now turn our attention to the Bg i~^i~'~f decay. The matrix element for the 
b si'^i^'-f decay can be obtained from that of the b si'^i~. For this purpose it is 
necessary to attach the photon to any charged internal as well as external line. When 
the photon is attached to internal lines, there will be a suppression factor of m^/m^ in 
the Wilson coefficient, since the resulting operators (dimension-8) are two order higher 
dimensionally than the usual ones (dimension-6). Thus the main contributions to the 
b s£+£^7 decay come from diagrams when photon is radiated from initial (structure 
dependent SD part) and final (internal Bremsstrahlung IB part) fermions. 
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In order to calculate the matrix element corresponding to SD part, the following matrix 
elements must be calculated in the first hand 

(7|s7^(l-75)6|5,) , 

These matrix elements can be parametrized in terms of two independent, gauge invariant, 
parity conserving and parity violating form factors [11, 16] as 



where and are the four-vector polarization and four-momentum of the photon, re- 
spectively, and q^ is the momentum transfer. Using Eqs. (3), (6) and (7), matrix element 
describing the SD part takes the following form 



MsD = Vt,vde^,^,e*^q>^k^\A,hH + A-,h>^^d 



+ i 



'e;{kq) - {e*q)k,] [b^I^^ + 52^7'^75^] } , 



(8) 



where 



Bi = Crf{q')-2Cr'^h{q'), 
When a photon is radiated from the final leptons, the corresponding matrix element is 



MiB = ^v,,y,:e/s,m,cr 



'10 



2pik 2p2k 



(9) 



Finally the total matrix element for the Bg 7 decay is equal to a sum of the Eqs. 

(8) and (9) 

M = MsD + MiB ■ 

The double differential decay width of the Bg — > £'^£~^ process, in the rest frame of the Bg 
meson, is 



dr 



M 



(10) 
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where and Ei are the energies of the photon and one of the final leptons, respectively, 
bar means that summation over spin of final particles is performed. The physical regions 
of E^ and Ei are determined from the following inequalities 



0< E^< 



TUB, - E^ 



E-y „ TUB, - 

2 - 2 



E^ Er« 



(11) 



where 



4m| 



is the lepton velocity. The term is infrared divergence free, the interference term 

2Re {M.sdM*ib) has integrable infrared singularity and only | Al/sl^ has infrared singularity 
due to the emission of soft photon. In the soft photon limit the Bs t^t~^ decay cannot 
be distinguished from the pure leptonic C.^l~ decay. Therefore, in order to obtain 

a finite result the Bs — > ^+£^7 and the pure leptonic Bg — > decay with radiative 

corrections must be considered together. It was shown in [17] that when both processes are 
considered together, all infrared singularities coming from the real photon emission and the 
virtual photon corrections are indeed canceled. However in the present work our point of 
view is slightly different, namely, instead of following the above mentioned procedure, we 
consider the Bg — > £+£"7 decay as a separate process but not as the 0{a) correction to the 
Bg — > decay. In other words, we consider the photon in the Bg — > t^l~^ decay as a 
hard and detectable one. Therefore, in order to obtain the decay width we impose a cut on 
the photon energy, which will correspond to the experimental cut imposed on the minimum 
energy for the detectable photon. We require the photon energy to be larger than 25 MeV, 
i.e., E^ > SmBj2, where 6 > 0.010. 

Using Eqs. (5) and (10) we get the following results for the Bg and Bg £+£"7 

decays, respectively: 



B{Bg i^r) = ^mljBjl \V,,V,t 



B{Bg ^ 

rl-4r 



rriz 



IT/ T/*|2 «^ _3 



'-'10 



(12) 



(27r) 



+ 
-4 



I ^' dxx\ [lAif + (1 - X + 2r) + (1^2!" + (1 - X - 4r) 

2rfB.Cl^ /'""'dxx2Re(^i)ln^ 
Jo 1 — V 



r-l-4r 



2v h 2H x] In 

X \ XX J 1 — V 



(13) 



where x = 2E^/mB, is the dimensionless photon energy and r = mj/m%^. 
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It follows from Eq. (13) that to be able to calculate the decay width, the explicit forms 
of the form factors g, /, gi and fi are needed. In further analysis, use will be made of 
the form factors that arc predicted by the light-cone QCD sum rules [11, 16], whose 
dependences, to a very good accuracy, can be written in the following dipole forms. 



1 GeV 



1 - 



2 ' 



(5.6 GeVf^ 
3.74 GeV^ 



1 - 



2 ' 



h<f) 



0.8 GeV 



1 - 



q^ 



2 ' 



(6.5 GeV)\ 
0.68 GeV^ 



40 



.5 GeV^ J 



1 - 



q^ 



30 GeV^ J 



(14) 



As is obvious from Eq. (2), to get numerical values of the branching ratios the value of the 
fourth generation CKM matrix element |Vj'fe\^r^| is needed. For this purpose we will use the 
experimentally measured values of the branching ratios B{B — > Xg'j) and B{B — > XcCUe)- 
In order to reduce the uncertainties coming from b quark mass, we consider the ratio 



R 



B{B ^ X,7) 



In leading logarithmic approximation this ratio can be written as 



(15) 



nf{mc)K{mc) \ Vc 



|2 ' 



(16) 



ch\ 



where rhc = rric/mb and the phase factor f{rh^ and C(q!s) QCD correction factor «;(mc 
[18] of 6 — c^v are given by 



f{rhc) = 1 - 8m^ + 8m^ - - 24m^ In(mc) , 



«(^c) = 1 - 
From Eqs. (16) and (15) we get 



31\ -.2\2 3 



TT 



4 ; 



(17) 



7r/(mc)fi:(mc) \Vcbf B{B Xs^) 



B{B ^ X,eu,) • 



(18) 



The model parameters can be constrained from the measured branching ratio of the 

B ^ Xgj decay 



B{B ^ X,7) = 



(3.21 ± 0.43 ± 0.27l|];}^) x 10"^ [19] 
(3.36 ± 0.53 ± 0.42 ± 0.54) x 10"^ [20] 
(3.11 ± 0.80 ± 0.72) X 10-^ [21] 
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In further numerical analysis, we will use the weighted average value B{B — > Xgj) — 
(3.23 ± 0.42) X 10'^ [22] for the branching ratio of the B X^i decay. 

Next constraint to the extended CKM matrix element comes from the unitarity condi- 
tion, i.e., 

^ = 1, 
^ = 1, 

Vu,v:, + k^k: + kok: + vt',vi = o . (i9) 

At this point we would hke to make the following remark: Charged-current three-level 
decays are well measured experimentally and they are not affected by new physics at leading 
order. For this reason, in the following discussions we will use Particle Data Group (PDG) 
constraints [10] for \Vud\-i |Ks|; iKftl and [Kb/Kbl- Using the weighted average for 

B{B X,7) and PDG constraint 0.38 < < 0.044 [10], from Eqs. (18) and (19) we 



1 '^us 




Ks 




IK 




iKsl 


1Kb 




[Kb 




IK 




iKfel 



obtain the following constraints 



0.011 < 



C7 KKs + CJ}''"Vt'iyt's\ < 0.015 , (20) 

0.03753 < IKK + KfeKsl < 0.043976 , (21) 

< |K|^ + IK,/ < 0.00492 , (22) 

0.998 < IKI^ + iKftl^ < 0.9985 . (23) 



3 Numerical analysis 

In this section we study the constraint on K'fcKs ^"<i using the resulting bound on [K'bKsl 
we investigate the influence of fourth generation on the branching ratios. Before performing 
this analysis we present the values of the main input parameters which we use in our calcula- 
tions: rub = 4.8 GeV, rric = 1.35 GeV, tub, = 5.369 GeV, tb^ = 1.64 x 10"^^ s. Furthermore 
in calculating the branching ratios the values of the Wilson coefficients Cf^(mfe), Cg^'^{mi)) 
and Cf^{mb) are needed. In leading logarithmic approximation C^^{mb) — —0.315 and 
Cf^{mb) = -4.642 (see [14, 15]). The analytic expression of C^^{mb) is given as 

Cg^irrib) = 4.227 +(7(m„s)(3Ci + C2 + 3C3 + C4 + 3C5 + C6) 

- ^(7(0, s){Cs + 3C4) - 1^(1, s){AC, + 4C4 + 3C5 + Ce) 

9 

where s = q^/rn^. The explicit form of the expressions for the functions g{mi,s) and the 
values of the individual Wilson coefficients can be found in [14, 15]. 

It should be noted that Bg — > £+£~7 can also receive long distance contributions which 
have their origin in the real intermediate states, i.e., J/'0, '0') ' ' "■ ^^e present work we 
neglect such long distance contributions. 

In order to study the influence of the fourth generation quarks to the rare Bg £^£^ 
and Bg — > £+£^7 decays we need the values of Vf^K ^.nd KfcKs- For this purpose we 
carry out the following analysis. We consider Vj^Ks ^t'bVfs be two independent 
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complex parameters. They are constrained by the unitarity condition (19) and the measured 
experimental value on branching ratio of the B Xg'-f decay (see Eq. (20)) which depends 
on mf. Thus for each value of there exists an allowed region in the (|VtbVj*| and |V^'6V^*s|) 
plane. 

It is not possible to solve the constraints Eqs. (20)- (23) analytically. Therefore, we chose 
a large number of random complex values for Vtb, Vts, Vfb, Vfs such that the constraints 
are all satisfied. For a sufficiently large number, the selected values would range over 
the whole solution space. In Figs. (l)-(3), to give an idea about the solution space, 
we present the allowed region in the plane of the absolute values of VtbV^l and VfbV^g 
at mt' — 200 GeV, rrit' — 300 GeV, and = 300 GeV, respectively. The upper left 
boundary of the allowed region is determined by Eq. (22) and (23) whereas the other 
three boundaries are determined by Eqs. (20) and (21). From these figures it can easily 
be seen that, when iVt/bV^jf^l ~ the value of |Vt(,l^*| is close to its SM value and is evenly 
distributed around 0.04. However, when the value of iVt/^V^T^I gets larger, up to a critical 
value around ~ 0.012, the allowed region for iT^ftV^*] becomes wider with the center 

fixed at 0.04. Beyond the critical point, the center of the region shifts toward smaller values 
with its width fixed. This behavior of the allowed range for continues until |V^t'6Vf*s| 

reaches a second critical value around |V^t'6V^t*s| ~ 0.04. When iV^'bVj?^! takes values beyond 
this value, the allowed region for iVt^Vj*! begins to shrink to zero with its center fixed 
around 0.012. Similar behavior is observed for all other values of mt'. After determining 
the physical regions for VtbV^l and Vt'bV^^g, the range of the region for the branching ratio 
of the Bs i^£^ and Bs £~^frj decays can be determined. Depicted in Figs. (4) 
and (5) are the dependence of the ranges of the branching ratios of the Bs — > A*"*"/^" and 
Bs — > T'^T~ decays on mt', respectively. It is evident from both figures that at lower values 
of mt' (200 GeV < mt' < 250 GeV) the branching value is distributed mainly around the 
SM value. This is to be expected since when m^ = mt, the contributions to the Wilson 
coefficients from the new physics effects are identical to those from the SM. In particular, 
after dividing Eq. (20) by Cf^ , Eq. (20) becomes a constraint on iVt^V^* + VfbV^J^l just as 
Eq. (21), and the branching ratio itself also depends on this combination only. Obviously, 
for larger values of mt' the possible departure from the SM prediction becomes substantial, 
signaling the existence of new physics beyond SM. 

It follows from these figures that there exists a region for VtbV*^ and VfbV^s where 
branching ratio attains smaller values than the SM prediction for all values of mf. This 
can be explained as follows. It follows from Eq. (12) that the branching ratio is proportional 
to 



B{Bs £+r) (X 



riSM 



1 + 



Vt'bVt 



t's '^10 



VtbVtl Cfn^ 



10 



\Vt,vz\ 



The ratio Cl'^"" 



I Cfo*^ is positive for all values of and increases with increasing mf 



From 



Eq. (24), we see that depending on whether |1 + A|, where 



A 



T/ T/* r^new 



is greater then, equal to, or less then one, the branching ration is greater then, equal to, 
or less then the SM value. The boundary separating these regions, |A + 1| = 1 can be 
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considered as a circle on the complex A plane with radius 1 and center at —1. Within this 
circle, the predicted branching ratio is smaller then the SM value, whereas out of this circle, 
it is bigger. 

The dependence of the branching ratio on mf for the Bg and Bg — > t"^t~7 

decays are shown in Figs. (6) and (7). Qualitative understanding of the dependence of the 
branching ratio on m^/ for the Bg — >■ t^l'^ {£ = /i, r) is analogous to the Bg {£ = 

^, r) case. The lower bound of x is chosen to be a; = 0.01 in order to get finite result for 
the B{{Bg — > £+£^7) decay (see Eq. (13)). It should be noted here that the branching ratio 
for the Bg fj,^fj,^j decay practically is insensitive to the variation of Xmin- For example 
when Xmin is varied from 10~^ to 10~^, the branching ratio changes only about ~ 5%. 

Our final remark is that fourth generation can give contribution to the rare B and 
K decays. For this reason constraints on fourth generation imposed by the p parameter, 
ArriK, £k, ^d{s) ~ ^d{s) Kiixing, — > n'^ui^, — >■ A*"*"/^" decays, must be considered 
together with h ^ decay. Early attempts in this direction were done in [22]. A more 
detailed analysis of constraints arising from rare B and K meson decays is presently under 
study. 

In conclusion, we have studied the effect of the fourth generation quarks to the rare 
Bg — > t^£~ and Bg — > £"'"£"7 decays. From the experimental result for the branching ratio 
of 5 — > Xg^ decay, the values of |Vt6V^*| and iV^'feV^*^! are constrained. We have observed 
that, depending on the relative phase between VtbV^g and VfbV^g, fourth generation quarks 
can enhance or suppress the SM result. As a result the branching ratios of Bg and 
Bg £+£~7 decays can be larger or smaller compared to that of the SM prediction. 
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Figure captions 

Fig. (1) The allowed regions for \VtbVZ\ and \Vt'bV;j at mf = 200 GeV. 
Fig. (2) The same as in Fig. (1), but at m*/ = 300 GeV. 
Fig. (3) The same as in Fig. (1), but at m*/ = 400 GeV. 

Fig. (4) The dependence of the branching ratio for the Bg — > /ji'^/jL" decay on the fourth 
up quark mass m^/. 

Fig. (5) The same as in Fig. (4), but for the Bs —>■ t^t~ decay. 
Fig. (6) The same as in Fig. (4), but for the Bg — > jJi^lJr^ decay. 
Fig. (7) The same as in Fig. (4), but for the Bg t^t~^ decay. 
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